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SIMPLE PROCESSORS OF STAR TRACKER COMMANDS FOR 


STABILIZING AN INERTIALLY ORIENTED SATELLITE 

By Robert D. Showman^ Brian F. Doolin^ 
and G . Michael Sullivan 

Ames Research Center 


SUMMARY 


The study develops processors which convert star tracker gimbal angle 
measijrements to satellite attitude control signals. An exact processor is 
nonlinear and complex because the gimbal angles vary as the inertial orienta- 
tion of the vehicle changes. In this paper ^ approximations to the exact pro- 
cessor are developed which greatly simplify the mechanization and yet allow 
satisfactory attitude control. 

The approximate processors are derived from a perturbation analysis in 
which only the first-order terms are considered significant . The first pro- 
cessor, called the ideal processor, is a direct mechanization of the first- 
order approximation to the exact equations and is the most complex to 
mechanize. It requires resolvers on the trackers* outer gimbal axes and the 
computation of the tangent and the sine and cosine of the difference of two 
angles. Two approximations to the ideal processor, called the partial pro- 
cessor and the constant processor, require less mechanization. The partial 
processor requires resolvers only on the trackers* outer gimbal axes. The 
constant processor only requires the mechanization of constants which are 
sign-controlled as a function of the outer gimbal angles. 

The use of the three processors to derive attitude control signals for an 
example satellite such as the Orbiting Astronomical Observatory is investi- 
gated. The functional requirements of the "coarse pointing mode" are imposed. 
It is shown that both the ideal and partial processor provide the required 
transient and steady-state performance over a commanded gimbal range that is 
limited only by the physical nature of the star tracker. Also, it is shown 
that resolver errors of 5^ in tbe mechanization of the partial processor have 
a negligible effect on the performance of the system. The constant processor 
is shown to provide a system with the required performance over a co mma nded 
gimbal angle range of ±6o° . 


INTRODUCTION 


An important class of satellites is one whose mission requires that the 
satellite orientation be fixed with respect to inertial space. Such a vehicle 
is usually termed "inertially oriented" or "inertially stabilized." Mission 
objectives for this type of satellite frequently require extremely accurate 
attitude control, thus a precision sensor is needed for determining attitude. 



One feasible sensor is the star tracker. The use of star trackers for 
attitude stabilization is the subject of this paper. 

The star tracker is an 
optical or light sensing device 
mounted in a gimbaled unit that 
has two degrees of freedom 
(fig. l) . The star tracker is 
locked onto a target star by 
sensing the misalinement between 
its optical axis and the line- 
of -sight (LOS) to the star and 
by using this error signal to 
drive the gimbals^ the star 
tracker controls its optical 
axis to be coincident with the 
LOS to the star. 

The attitude of a satellite 
can be related^ in general^ to 
a reference frame formed by the 
directions to two stars . If 
two star trackers are pointed 
at known guide stars and if the 
orientations of the star trackers with respect to the vehicle are known^ the 
attitude of the satellite can be determined. Conversely, if the vehicle is in 
a desired attitijde, the relative orientations of the two star trackers with 
respect to the vehicle can be calculated in terms of their gimbal angles. If 
the calculated gimbal angles are the commanded angles, the deviation of the 
gimbals from their commanded angles then provides a measure of the satellite's 
deviation from its prescribed attitude. At least three of the four gimbal 
angle errors from the two star trackers are required to determine the vehicle 
attitude error. Since the gimbal angle errors are generated about axes which 
are, in general, neither alined to the control axes of the vehicle nor mutu- 
ally orthogonal, they must be processed to be useful as control signals. 

Computation of the satellite attitude errors from the gimbal angle errors 
is quite complex if the exact equations are used. An approximate relation- 
ship is desired that simplifies the computation. If, however, the satellite 
attitude and guide stars for the trackers are to be arbitrary, then the gimbal 
error processor must provide useful control signals over a wide range of star 
tracker gimbal angles. Yet, if the processor is to be used on board a 
satellite, it should be as simple as possible. 

The use of star trackers to derive control signals for satellite stabili- 
zation or attitude control has been studied in references 1 to 3. Meyer 
developed a scheme that can be used either to stabilize the attitude of the 
satellite or to control it through large-angle reorientation maneuvers. Since 
the method is general, computation of the control signals is complex and 
requires a computer on board the satellite. The processor developed in ref- 
erences 2 and 3 is much simpler to mechanize and does not require a computer. 
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It is^ however^ applicable only for attitude stabilization. The scheme is 
derived on the basis of a first-order approximation of the satellite from its 
nominal or prescribed attitude. Each star tracker provides an independent 
estimate of the satellite’s error in attitude. Since a minimum of two 
trackers is required^ the individual estimates are averaged to provide an 
approximation to the attitude error. Although its mechanization is simple^ 
the scheme does not allow attitude stabilization of the satellite over the 
complete range of gimbal angles^ since instabilities can occur when the com- 
manded gimbal angles become large. 

This paper considers three related processors of gimbal angle errors . 

They are like the processor in references 2 and 3 ^ in that they come from a 
perturbation analysis where only first-order terms are considered significant. 
They differ^ however^ in that they process gimbal angle errors from pairs of 
star trackers rather than from each tracker independently. Of the three 
methods of processing the gimbal angle errors, the first, called the ideal 
processor, is a direct mechanization of the first-order approximation to the 
exact eq-uations. The other two schemes, called the partial processor and the 
constant processor, are variations of the ideal processor that simplify the 
mechanization and eliminate the need for a computer. The ideal processor 
does require more computation than the simplified forms. The characteristics 
of the system using each of the processors are developed so that the choice of 
a particular processor for a specific application can be made by the system 
designer . 

The first section of the report discusses the problem of determining a 
satellite’s attitude by means of star trackers. The pertinent attitude per- 
turbation equations are derived and various conditions that simplify their 
form are examined. The three processors are then defined and their peciiliar- 
ities discussed. The next section concerns the design and performance analy- 
sis of a system using each of the three processors. Here the functional 
requirements and basic stabilization system parameters of the Orbiting Astro- 
nomical Observatory are used for illustration. Finally, a section is devoted 
to the results of a computer simulation of the example stabilization system 
using each of the three processors . 
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ATTITUDE DETERMINATIOW VIA STAR TRACKERS 


The role of the star trackers in the control system is depicted in 
sketch (a). The star trackers, represented by one block in the sketch, are 
considered to be locked onto selected target stars. When the vehicle is in 
its desired attitude, the trackers point away from the satellite in the direc- 
tions corresponding to the gimbal angles which are selected on the ground and 
stored in the vehicle as commands. But since various disturbances move the 
vehicle while the trackers still point at their targets, the gimbal angles 
differ from their commanded values. The sketch illustrates that these angle 
errors are detected and processed to form an estimate of the vehicle's devia- 
tion from its desired attitude. Finally, the sketch shows that this estimate 
is used by the controller to activate a motor or jet, which, in turn, drives 
the vehicle back to its desired attitude. 
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Sketch (a) 


The forms of the geometry matrix, N, and the processor matrix depend on 
the vay the star trackers are mounted on the vehicle. Figure 2 illustrates 
the mounting arrangements considered in this report, defining them relative 
to the vehicle's control axes (l^ - roll, 2^ - pitch, 3^ - yaw). The trackers 
are shown with their gimbal axes in the null position and their optical, 
inner, and outer axes labeled 1, 2, and 3, respectively. The superscript indi- 
cates the tracker number; thus, 1^ is the optical axis of tracker number 2. 
Trackers 1, 2, and 3 are mounted to the vehicle similarly: their outer gimbal 

axes ( 3 ^) are alined to the roll axis ( 1 ^) of the vehicle. Tracker 4 is 
mounted differently: its outer gimbal axis, 3^, is parallel to the pitch axis 

(2^) of the vehicle. 

The orientation of a star tracker with respect to the vehicle is speci- 
fied by two gimbal angles. The angles for tracker number 1 are defined in the 
drawing on the right in figure 2. The inner and outer gimbal angles are 
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Figure 2.- Illustration of gimbai angles and arrangements of star trackers on the spacecraft. 


designated^ and Jxy respectively. To define a coordinate system for the 
star tracker, a third angle (o-i) is introduced which refers to angular rota- 
tion about the optical axis. The angle aj_ cannot be measured by the star 
tracker. Angles for the other trackers are similarly defined. 

The motion of the satellite can be expressed both in the coordinate sys- 
tem defined in the vehicle and in the coordinate system defined in the tracker. 
A simple way of describing the motion is to express the angular velocity of 
the vehicle ^av^ ^ 3 v) terms of the gimbai angle rates of the star 

tracker (ai. Pi, 7i)« In order to equate the gimbai angle rates to the vehi- 
cle angular rates, it must be assumed that the sensor tracks its star per- 
fectly. The sensor's error signal, which this assumption ignores, can easily 
be added later if necessary. The transformation between tracker gimbai rates 
and vehicle angular rates is derived in appendix A. The rate transformation, 

R, for each of the trackers indicated in figure 2 is given in table I. Since 
the control system is to be used for stabilization and the motion about any 
nominal attitude can be considered small, a perturbation analysis is appropri- 
ate. Therefore, the vehicle angular rates ^^y^ ^sv^ approximated by 
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the Euler angle rates (cp^ Gy * 1 ^) . Since the motion is assumed to be small^ 
trigonometric elements in the transformation R can be considered constants 
that are evaluated at the nominal or commanded angles. The transformation R 
then provides a first-order approximation to the relationship between the 
deviation of the gimbal angles from their commanded values and the deviation 
of the satellite from its prescribed attitude. For example^ the first-order 
approximation for tracker number 1 is 




0 C7i/cpi -S7i/cpi 



Api 

= 

0 S7i C7i 


A0 

A7i 


1 -C7itpi S7itpi 


i Ailr 


The attitude error of the vehicle could be determined from equation (l) 
if each of the error signals Acx^i^ and A71 were measured by the star 

tracker. Since the star tracker cannot sense angular motion about its optical 
axis^ Ao-i is not measurable. A gimbal angle error signal from a second star 
tracker must then be substituted for Z^i to determine the vehicle attitude. 
The third measurable signal could be either an outer gimbal error signal (A7) 
or an inner gimbal error signal (AP) from any of the remaining trackers. For 
example^ if the inner gimbal error from star tracker 3 (APs) is substituted 
for A:x-i^ the following transformation is obtained: 


APi 

A 7 i 

AP3 


0 S7i 

1 -c7itPi 

0 -C7s 


— 


— 

C7i 


Acp 

S7itpi 


A0 

S73 


At 


( 2 ) 


The attitude of the vehicle can now be estimated from the gimbal angle 
measurements (if they form an independent set) by inverting equation (2). 

Equation ( 2 ) gives a form that matrix N can take. Several forms aris- 
ing from different combinations of gimbal error signals from two trackers are 
given in table II; many other combinations are also possible. For example^ 
one error signal from each of three star trackers could be used to estimate 
the error in vehicle attitude. However^ the use of gimbal errors from three 
trackers has a distinct disadvantage when reliability is considered. Conse- 
quently^ gimbal error signals from only two trackers are considered. The 
giinbal error combinations listed in table II are representative of ail the 
possible cases. 
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If the gimbal freedom of the star trackers Is restricted to small angles 
and if the trackers are properly mounted to the vehicle, the gimbal error sig- 
nals directly provide an accurate estimate of the attitude error. However, 
it is usually desirable to be able to command large gimbal angles to track the 
available guide stars. Therefore, the gimbal axes are, in general, not alined 
to the control axes nor do they fonn an orthogonal set. A transformation is 
then required to convert the gimbal error signals into an estimate of the 
vehicle attitude error. Three classes of such transformations, or processors, 
are defined. They are referred to as (l) the ideal processor, (2) the partial 
processor, and ( 3 ) the constant processor. 


Ideal Processor 

The ideal processor is defined as the first-order approximation relating 
the vehicle control signals to the gimbal angle errors, For tracker pair 1,3 
when the gimbal errors APi, A/i, and A ^3 are used, the processor is the 
inverse of equation (2) and is given by the following equation; 


€9 


s(7s-7i)tPi c(7i-7s) 


1 

0(71-73) 


S73 

073 


0 

0 


— ' 


, 

-tpl 


APi 

-C7i 


A7i 

S7i 


Aps 

_ 





( 3 ) 


The ideal processor for the other combinations of gimbal angle errors is given 
as the M matrix in table II. Each is the inverse of the corresponding N 
matrix. 


The processor (M) could be mechanized on the satellite in at least two 
ways. The first method would be to mechanize each of the matrix elements with 
a potentiometer whose values would be programmed for each desired attitude. 

The potentiometers could be programmed either manually from a ground station 
or automatically on board the satellite. In either case, additional telemetry 
channels are required. Manual programming from a ground station has the 
advantage that an on board memory unit is not required but the disadvantage 
that changes in attitude can be made only when the satellite is over a ground 
station. Automatic programming of the potentiometers has the advantage that 
changes in the desired attitude can be made at any position in the orbit but 
the disadvantage that a memory unit is required. The size of the memory unit 
depends on the number of changes in desired attitude between ground stations 
and the number of tracker pairs that must be programmed. The second method 
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would be to mechanize the trigonometric elements actively. Although such 
mechanization is more complex than that with potentiometers^ a storage unit 
and additional telemetry channels are not required. Also^ attitude can be 
changed at any position in orbit. The theoretical investigation that follows 
is applicable to either mechanization. 

The complexity of the ideal processor varies with the particular gimbal 
error combination selected. For example^ the ideal processor for tracker 
pair 1^2 - if and AP^ are used - is 



and has the same form as equation (3). The processors in equations (3) and 
( 4 ) are clearly simpler than the other ideal processors listed in table II. 

The simpler processors have not only two zero elements but also simpler 
trigonometric functions for the nonzero elements. 

The complexity of the ideal processor depends on both the relative aline- 
ment of the paired star trackers and on the combination of gimbal errors 
chosen to estimate the attitude error. Although many combinations were 
studied^ only representative cases are discussed. The variation in complexity 
due to the relative alinement of the paired trackers is revealed by comparing 
the processors for tracker pairs 1^3 and 1 ^ 4 . Since both use a gimbal error 
combination of two inner and one outer ^ the simpler processor of tracker pair 
1^3 clearly indicates the preferred alinement occurs when the outer gimbal 
axes are parallel. The variation in complexity due to the combination of gim- 
bal errors then is exhibited by comparing the two processors for tracker pair 
1 ^ 2 . The simpler processor is obtained when a gimbal error combination of 
two inner (APi^ AP^) Q'^d one outer (A71) is used. In summary^ therefore^ the 
simplest ideal processor is obtained when the paired star trackers are mounted 
to the vehicle so that their outer gimbal axes are parallel and when a combi- 
nation of two inner and one outer gimbal error signals is chosen. Mechanizing 
this processor insures that^ to a first order^ the control channels are 
decoupled. 
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Partial Processor 


Mechanization of even the simplest form of the ideal processor for use 
on hoard a spacecraft is relatively complex due to the appearance of such 
functions as the sine and cosine of the difference of angles and the tangent. 
It is usually desirable to obtain a simpler form that still provides the nec- 
essary performance. Simpler forms are obtained by considering approximations 
to the ideal processor. If^ to reduce complexity^ only simple sine and cosine 
terms in the ideal processor are mechanized by using resolvers, while the more 
complex elements are made constant, the transformation in equation (3) can be 
written 


Mp = 


Vxi 

1 

VX3 

<3-X3S7s 

0 

-di3C7x 

(3.13C73 

0 

<ii3S7i 


( 5 ) 


where Vix is the constant for -[s(7i-73)/c(7i-73) ]tPi, V13 is the constant 
for [ -tPx/c(7i-73) ] , and dx3 is the constant for [ l/c(7x“73) ] • The con- 
stants (vxx^ dxs) are independent of the gimbal angles. Eq.uation ( 5 ) is 

designated the partial processor. The product of the partial processor and 
the geometry matrix is 


1 (^XiSri-C/itpi-VisCrs) (viiC7i+S7itPi+Vi3S73) 


MpN = 


0 


d.xsc(7i-73) 


0 


0 


0 

'ii3c(7i-73) 


( 6 ) 


It follows from equation (6) that this processor insures the independence of 
only two control signals and e,) rather than all three as obtained with 

the ideal processor. However, the^loss of independence of is a small 

price to pay for the simpler mechanization. Equation (6) also shows that the 
gains of £q and through the partial processor are functions of the com- 
manded outer gimbal angles, whereas their gains through the ideal processor 
are constant. Consequently, the dynamic response of the system using the 
partial processor varies as a function of the commanded outer gimbal angles. 


The partial processor is the processor that provides two independent 
control signals when only the sine and cosine terms of the ideal processor are 
mechanized. It is derived from only the simplest ideal processor (table II). 
To use this processor, the inner gimbal error signal from one tracker is 
passed through a resolver on the outer gimbal of the second tracker of the 
pair; thus, APx passed through the outer gimbal resolver of 73, and APs 

is passed through the outer gimbal resolver of 7x* 
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Constant Processor 


Although the partial processor allows simple mechanization^ another 
approximation to the ideal processor allows even simpler mechanization. This 
form^ designated as the constant processor^ is obtained by letting each of the 
nonzero elements of the ideal processor be constant. The elements of the 
matrix are not parameters to be evaluated at each commanded set of gimbal 
angles. They are to remain constant at least in magnitude over the whole 
range of allowed values of the angles. However^ the matrix elements have dif- 
ferent signs in different regions of the range of angles. These sign changes 
will be discussed later. The constant processor for tracker pair 1^3 


^11 

1 

^13 

rsi 

0 

^23 

rsi 

0 

^33 


( 7 ) 


where r^x is the constant for -[ s(7x-73)/c(7x-7s) ] t3x^ is the constant 

for s73/c(7x-7s) ^ etc. Only seven of the nine matrix elements must be 
mechanized. The product matrix is 


MeN = 


1 (riiS7i-c7it(3i-nsC73) 

0 (r2iS7i-r23C7s) 

0 (r3lS7i-r33C73) 


(3^‘iiC7i+s7i'bPi+ri3S73) 

(r2iC7i+r23S73) 

(i:’3lC7l+r33S73) 


(8) 


and shows that the pitch and yaw error signals are independent of the motion 
about the roll axis. 

The constant processors derived from the simplest ideal processor 
(table II) not only simplifies the mechanization but also simplifies the anal- 
ysis. The mechanization is simplified because only seven of the nine elements 
must be considered. The constant processors which might be derived from ideal 
processors other than the simplest require that all nine elements be mecha- 
nized. The analysis is simplified because the control signals are partially 
decoupled (see eq. (8)). 


Indeterminant Condition 

Since the gimbal axes vary as a function of the relative orientation of 
the star tracker and vehicle^ the error signals do not^ in general^ form an 
orthogonal set of measurements and can^ in fact^ form a dependent set of mea- 
surements. When this occurs^ the attitude of the vehicle is unobservable. 
This situation is termed the indeterminant condition. 
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Figure 3. - Illustration shoving when indeterminant condition would occur. 


The indeterminant condition exists when the determinant of the N matrix 
vanishes. The indeterminancy is therefore independent of the processor. From 
the determinants given in table II (D12, Di3^ can be seen that 

the condition is simplest to ascertain when the star trackers are used in a 
manner that provides the simplest ideal processor. For example^ the indeter- 
minant condition occurs for tracker pair 1^3 when D13 = c(7i-73) = 0^ for 
tracker pair 1^2 when ^ ^(71-73) = 0^ and for tracker pair 1^4 when 

Di 4 = s7iS74-c74tPx = 0 . The geometrical interpretation of the indeterminant 
condition for the preferred design (D12 and D13) is that the optical axes of 
the paired star trackers and the outer gimbal axes are coplanar. The situa- 
tion would occur if the direction associated with point C in figure 3 were 
contained in the plane defined by the vehicle and the targets stars A and B. 
This interpretation provides a simple criterion for the selection of stars by 
ground operations to avoid the condition. 

The vehicle cannot be operated within the neighborhood of the indetermi- 
nant condition because the error signals become less and less sensitive to 
some motions. For example^ the indeterminant condition for tracker pair 1^3 
occurs when 0(71-73) = 0 or^ equivalently^ when I71-73I = 90 °. Therefore^ 
the region in which vehicle operation is prohibited is defined as 

90° - ^ l7i-7sl ^ 90° + (9) 

where defines the magnitude of the restricted region. A graphical 

interpretation of the restricted region is indicated by the shaded area in 
sketch (b). 







Sketch (b) 

Simlariy^ the indeterminant condition for tracker pair 1^2 occurs -when 
3(71-72) = 0 or )7i“72| = 0^. The restricted region is 


1/1-72 1 ^ A7j^ 

and is indicated hy the shaded area in sketch (c). 
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ANALYSIS AND DESIGN 


The use of the three processors in deriving control signals for attitude 
stabilization of an example satellite is investigated in this section. An 
example of a satellite that might use the processors developed in the previous 
section is the Orbiting Astronomical Observatory (OAO) whose mission objec- 
tives require that the experiments on the satellite be directed at arbitrary 
points in the celestial sphere. When the vehicle is pointed at an arbitrary 
target^ the attitude is defined by gimbaled star trackers and is controlled 
by three mutually orthogonal motor inertia -wheel combinations. Both the sta- 
bility and performance of the example system using each of the processors to 
derive control signals will be investigated. Although the example satellite 
uses momentum transfer for control^ the use of the processors is not 
restricted to such a system. 
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Description of Control Systems 



Figure 4.- Attitude control system using star trackers I and 3 . 


Figure 4 is a block diagram of the control system- The gimbal error 
signals from tracker pair 1^3 shown being processed to provide an estimate 

of the attitude errors. The processor indicated on the diagram could be any 
one of the three discussed in the previous section^ ideal^ partial^ or con- 
stant- The equations provided in the diagram assume the following conditions: 
(l) The mass distribution of the vehicle is such that the moments of inertia 
about any three orthogonal axes are equal; ( 2 ) the linearized equations val- 
idly describe the motor inertia-wheel combination; ( 3 ) the motor torque is 
limited; and (4) gyroscopic coupling due to inertial wheel rotation is negli- 
gible. A lead compensator is introduced because without it the system would 
be highly oscillatory. 

The motor nonlinearity shown in figure 4 is represented by a describing 
function (ref. 4). This representation is valid since the motor-vehicle com- 
bination behaves as a low-pass filter. The describing function is a pure gain 
which decreases monotonically from its maximum value of unity - obtained with 
the input in the linear region - to zero as the input increases. 


Ideal Processor 

Stability . - Mechanization of the ideal processor insures that^ to a first 
order ^ the error signals^ and are independent; that is^ each error 

signal is only a fionction of the attitude error about its corresponding axis. 
For example, the product of the ideal processor M and the geometry N is 
the unity matrix. The control channels are therefore uncoupled and system 
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stability is determined by investigating each control loop independently. 

That this is true can be seen by deriving the characteristic eq.uation (CE). 

CE = DET(I + GMN) = 0 (ll) 


■where 
MN = I 


( 12 ) 


G = 


Si 

0 

0 


0 

§2 

0 


0 

0 

gal 


(13) 


g,- = 


NjLK(TiS + 1) 


i = 1. 2, 3 


"“i s(t2S + 1)(%S + 1) 

= describing function for the ith channel 


(U) 


K = 


(15) 


Equation (ll) can be written as follows: 

CE = (l+ gi)(l+ g2)(l+ ga) 


= 0 


(16) 


Since the motion is damped if all roots of the characteristic equation have 
negative real parts^ it is necessary to show that the roots of each factor in 
equation (l6) have negative real parts. A general expression for the factors 


IS 


(l + gj) - 0 i - 1^ 2^ 3 


(17) 


Substituting equation (l4) into equation (l?) gives: 


T2 Tr, ' 




(18) 


where 


= 10T2 


l6 




The parameter plane method (ref. 5) is convenient for determining sta- 
bility. This method allows the system’s stability to be expressed as a func- 
tion of two variable parameters. If the two variable parameters are defined 
by the following equations^ 


B1 - 1/t2 
B2 = R^K 


(19) 


the characteristic equation can be written as follows: 


s3 + 


Bin- 

T- 


S^ + 


B1 10B2 


m 


'm 


s + 


B1B2 


'm 


( 20 ) 


By substituting ® equating the real and imagi- 

nary parts of the resulting expression for equation (20) to zero^ the variable 
parameters can be expressed as functions of damping (^) and natural frequency 
(tOj^) for a particular value of Thus^ contours of constant ^ can be 

mapped as functions of B1 and B2 with the ^ = 0 contour defining the 
stability boundary. 





Sketch (d) 


The stability regions on the parameter 
plane are obtained by determining both the 
real root boundaries (a = 0 and a = oo) and 
the complex root boundaries (^ = 0). Real 
root boundaries (ref. 5) occur at Bl = 0 
and B2 = 0 and are plotted along with the 
complex root boundaries in sketch (d). If 
the boundaries are properly shaded^ ^ the 
number of stable roots (roots with negative 
real parts) can be determined for any combi- 
nation of Bl and B2. The numbers in the 
sketch indicate the number of stable roots 
in each region. Since the characteristic 
equation is third order^ the sketch shows 
that the system is stable, provided both 
Bl and B2 are positive. An expanded plot 
of this first quadrant is given in figure 5* 


The stability of the system using the ideal processor is determined by 
examining the range of the variable parameters (Bl, B2). It is observed from 
equation (l9) that Bl and B2 are always positive because and K 

are always positive. Since the describing function % varies from zero to 
one, the operating point of the system varies along a constant Bl line from 
B2 0 (large errors) to B2 = K (linear region) but is always in the first 
quadrant. The system is, therefore, stable. 


^The shading of stability boundaries is opposite to that suggested in 
reference 5- 


17 




Figure 5.- Parameter plane for ideal and partial processor. 


Practical considerations require that certain gimbal angles should not 
be commanded. The analysis just concluded Implies that the system is stable 
for any range of commanded gimbal angles. Although this is true, mechaniza- 
tion of the ideal processor requires the term 1/0(71-73) to be mechanized. 
Since [1/0(71-73)] “ as l7i-73t 90 °, an upper limit must be imposed on 

the gain term. This limit has already been established in the discussion on 
the indeterminant condition and was defined as the region of restricted 
operation in the 71, 73 plane. 

The simultaneous use of multiple tracker pairs in the control loop may 
be desirable for reliability pirrposes. The advantage of using more than one 
pair of trackers is that the malfunction of one tracker will not result in 
the loss of vehicle control. If several pairs are used simultaneously, the 
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estimated attitude error signals from each pair of trackers may he averaged. 
For example^ if the loop gains for tracker pair 1^2 and 1^3 are^ respectively, 
B2 i 2 B2 i 3, the loop gain for the multiple pair system is 


B2 


av 


B2i2 + B2i3 
2 


( 21 ) 


Since the B2 for any individual tracker pair is positive, the average is 
positive. The system using multiple pairs of trackers is therefore stable. 

Performance . - The system using the ideal processor has been shown to be 
stable. It is now necessary to choose the parameters B1 and B2 that will 
provide the desired transient and steady-state performance - 

The OAO telescope is alined with the roil axis of the vehicle. Since the 
telescope must point accurately in some specific direction, two stringent 
steady-state conditions are imposed on the pitch and yaw channels. The first 
condition, called Pointing Accuracy, requires that the roll axis of the vehi- 
cle be pointed to within ±1 min of arc. The second condition, called Drift 
Rate Accuracy, requires that the desired attitude be held to within 15 sec of 
arc for 50 minutes of time. 

The steady-state errors are obtained by applying the final value theorem 
to the pitch and yaw transfer functions when the torque input is a step func- 
tion. The pointing error and drift rate error are 



where and (aVso are initial wheel speeds for pitch and yaw motors; "ss" 

refers to steady state; t^ and t3 are external torques about the pitch and 
yaw axes; and PE and DR are the maximum acceptable pointing error and the 
maximum acceptable drift rate error, respectively. In addition, it is assumed 
that the steady-state motion is in the linear region so that = 1. 

Equations (22) and (23) rewritten as 


^2]pE ^ 




CO, 


W20 




V30 


2-1 1/2 


( 24 ) 


B2 


DR 



(25) 
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External torque ( |t2l = |t j|), dyne-cm 


Figure 6.- Steady-state gain requirements - ideal 
and partial processor. 


specify the ininimuin gain K (or 
minimum B2) necessary to satisfy 
each of the steady-state condi- 
tions. The minimum gain (K) cal- 
culated from equations (24) and 
(25) is plotted in figure 6 for 
the conditions listed in 
table III. For example^ if 
l^w 2 ol ” percent of the maxi- 
mum wheel speed and if 
|t 2 | = 1000 dyne-cm^ the minimum 
gains (K) necessary to satisfy 
the pointing accuracy and the 
drift rate accuracy are^ respec- 
tively^ K = 2.7 and K = 3*0* 
Therefore^ the minimum gain that 
satisfies the drift rate specifi- 
cation also satisfies the point- 
ing accuracy specification. This 
minimum gain is equivalent to the 
minimum value of B2 in the lin- 
ear region (Nj_ = 1) . A value of 
B1 must then be chosen in the 
acceptable range of B2 so that 
the vehicle exhibits the desired 
transient performance in the 
pitch and yaw channels; that is^ 
a B1 and B2 (n^_ = l) are 
selected from the parameter plane 
in figure 5 so that the desired 
damping ratio and natural 
frequency are obtained. 


Since the observatory is alined to the roll axis of the vehicle^ perfor- 
mance in roll is not so critical as in pitch and yaw. However^ the procedure 
for designing the roll control channel is identical to the procedure just 
described for the pitch and yaw channels. 


Partial Processor With Perfect Mechanization 

Stability. - The characteristic equation (CE) for the system using the 
partial processor is derived as follows: 

CE = DET(l + GMpN) = 0 (26) 

The expanded form of the characteristic equation, with tracker pair 1,3^ is 

(1 + gi)[l + g 2 <i 132 ( 71 - 73 ) ][1 + gs<ii3c(7i-73) ] = 0 (27) 
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where the are given by equation (l4). Since the stability of the system 
is determined by investigating the roots of the individual products, the fol- 
lowing generalized form is defined. 


(1 + « 0 


where 


( 28 ) 


K 


SI 


1 i = 1 

d-i3c(7i-73) 1 = 2, 3 


(29) 


Equation (l4) can be used to write equation (28) as 


Tg T, 


1 . 1 ^ + 10WlKslK ^) ^ ^ KNjKsi ^ ^ 

T T T 

2 m m 




(30) 


where = 10 x 2 * 

The parameter plane is used to determine stability. If the variable 
parameters are defined as follows, 


B1 



B2 = 


(31) 


equation ( 30 ) becomes identical to equation (20). Therefore, the regions of 
stability defined on the parameter plane for the partial processor are identi- 
cal to the regions defined for the ideal processor (sketch (d) and fig. 5 )- 
The difference between the behavior of the system using the processors lies in 
the variability of the parameter (cf. eqs, ( 3 I) and ( 19 )). 

The system using the partial processor is then stable if the parameters 
B1 and B2 are always positive. Clearly, the parameter B1 is positive. 

Since and K are positive, the condition B2 > 0 is satisfied if Kgj_ > 0. 
The condition is satisfied for the roll channel” (Kg^ = 1) - For the pitch and 
yaw channels, the condition > 0 requires that the sign of dis be 

controlled as follows: 


di3 ^ 0 if |7i-73lS90° ( 32 ) 

The parameter B2 is then positive for ail three control channels and the 
system is stable. 

The following observations are made from the preceding analysis. First, 
the two constants Vix and V 13 in the partial processor for tracker pair 1,3 
do not affect the stability of the system. They do, however, influence the 
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amount of pitch and yaw motion that is coupled into the roll channel (see 
eq. (6)). To give Vn and V13 values other than zero in the hope of mini- 
mizing the coupling leads either to complex logic or to computation equiva- 
lent to mechanizing the ideal processor. Henceforth^ they will be set equal 
to zero. Second^ the technique of coiribining outputs from multiple pairs of 
trackers when using the ideal processor may be applied successfully when using 
the partial processor. 


Performance. - The parameters must now be chosen so that the system will 
exhibit the desired transient and steady-state performance. Using the defini- 
tions of steady-state error from the previous section^ inequalities analogous 
to equations ( 24 ) and (25) are 


B 2 

1 

= 

Lx. 


(ho - 


'PE 

S2 

IPE 

PE 

Lv ° 

B 2 

1 

= 


s — 



'dr 

S2 1 

DR 

DR 

.vv 





1/2 


(33) 

(3^) 


where 

^S2 "" ^S3 ^^9)) 


Equations (33) and ( 34 ) provide the minimum loop gain necessary to 

satisfy, respectively, the pointing accuracy requirement and the drift rate 
requirement. The system will satisfy both steady-state requirements if the 
magnitude of the gain is chosen at least as large as the larger of the gains 
specified by equations (33) and (3^)^ that is. 


B2 ^ max 


KK. 







( 35 ) 


where 


Wj_ = 1 (linear region) 


The transient performance of the system varies with the relative orienta- 
tion of the star trackers and vehicle because of the variation in loop gain 
(B 2 ). The desired transient perfozmiance must he achieved over the complete 
range of parameters. Assuming that the minimum B 2 is chosen in the range 
specified hy equation (35)^ equation (3I) becomes 


Kdi3 


®^min 

[c(7i-73)]min 


(36) 


where the minimum value of 0(71-73) is dictated by the magnitude of the 
restricted region. The maximum value of B 2 is 
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®2jnax - ^^^13 - 


B2. 


min 


[c(7i-73)]mi 


min 


(37) 


A value of B1 is chosen from the parameter plane (fig. 5) so that acceptable 
transient performance in pitch and yaw is achieved over the complete range of 
B2, The pitch and yaw channels then exhibit both the desired transient and 
steady-state performance. The desired performance in roll is not so stringent 
as pitch and yaw and is easily achieved. 


The system using the partial processor for other tracker pairs exhibits 
performance identical to that just found for tracker pair 1,3« Also^ the 
design procedure is identical. The only difference between tracker pairs 
occurs in the definition of Kgj_. If the paired trackers are mounted on per- 
pendicular faces of the vehicle^ the gain == paired 

trackers are mounted on parallel faces of the vehicle^ ihe gain 

^si ” ^jk^(7 j -7k) • 


Partial Processor With Imperfect Mechanization 


General . - The partial processor^ if perfectly mechanized^ provides a sys- 
tem that stabilizes the satellite over a range of commanded gimbal angles 
limited only by the possible angular rotation of the gimbals. Since errors 
in mechanization exists it is desired to determine their effect on the stabil- 
ity and performance of the system. Mechanization of the partial processor 
involves passing the gimbal error signals through resolvers and through fixed 
gain elements (pots or amplifiers). The fixed gain elements are assumed to be 
perfectly mechanized and the errors in mechanization exist only in the 
resolvers. Imperfect mechanization results from either the misalinement of 
the resolvers or an inaccurate readout of the gimbal position and would appear 
as a resolver angle of 7 + A 7 where A 7 represents the error. The errors 
are assumed to be less than 5 ^ so that only the first-order terms are 
significant. 

The partial processor for tracker pair 1^3 with errors in mechanization 
is 


Mpg = di3s(73+A7s) 

di3c(73+A73) 


1 

0 

0 


-di3c(7i+A7i) 

di3s(7i+A7i) 


( 38 ) 


The product of the processor matrix (Mpg) and the geometry matrix (W) would 
show that the errors in mechanization have the effect of coupling the pitch 
and yaw channels. The effect of the errors on the stability and performance 
of the system is now investigated. 
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stability . - The characteristic equation (CE) for the system using the 
partial processor for tracker pair 1,3 with errors in mechanization is 


CE = DET(I + GMpjjW) =0 (39) 

The expanded form of equation (39) is derived in appendix B (B9) and is 

CE = (1 + gNi)(l + B^g + BpgS) = 0 (40) 


where 

g = (41) 

s(t2S + 1)(%S + 1) 

Bp = (N 2 + N3)di3c(7i-7s) + (ii3(c73S7iN2 - s73C7iNs)5 (42) 

Bj) = W2N3di3c(7i-73)[di3c(7i-7s) + di3s(7i-7s)5] (43) 

6 = (A 73 -A 71 ) (44) 


Parameters Bp and Bp are also derived in appendix B ((B14) and (BI 6 ) ) and 
are obtained by considering only the first-order error terms as significant. 
The parameter 5, defined by equation (44), is a measure of the error in mech- 
anization. The coupling introduced by the mechanization errors precludes 
factoring the second -order polynomial in g (see eq. (40)) into two first- 
order polynomials. (Compare equation (27) for perfect mechanization.) 


The stability of the system is determined by investigating the factors 
in equation (40). Since the roots of the factor (l + gNi) = (l + gi) have 
already been shown to be stable, it is only necessary to Investigate the 
polynomial 



1 + Bpg + BpgS = 0 (45) 

System stability can be expressed as a func- 
tion of the two parameters and by 

using the parameter plane method. The real 
and complex boundaries as well as the number 
of stable roots for each region are shown in 
sketch (e). Since the polynomial of equa- 
tion (45) is sixth order when expanded^ any 
combination of Bip and Bj^ that remains in 
the first quadrant but below the ^ = 0 
curve provides a stable system. An expanded 
plot of the first quadrant is given in 
figure 7* 
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Bt 

Figure l.- Parameter plane - partial processor - 
trackers 1^3. 


The stability of the system 

is determined by calculating the 

range of the parameters B,p and 

Bq. If the actual parameters are 

contained within the stable 

region indicated in sketch (e)^ 

the system is stable. Both the 

minimum and the maximum Bj^ must 

be calculated for a given Bip. 

If Bn is less than the 

Bmax 

^ = 0 boundary and if Bpj . >0^ 

the system is stable. 

For comparison purposes^ the 
stability of the system with per- 
fect mechanization is determined 
before imperfections are intro- 
duced. The system has already 
been shown to be stable. The dif- 
ference between the previous anal- 
ysis and the present analysis is 
in the definition of the equations 
used to express stability. For 
perfect mechanization^ equa- 
tions (42) and (43) become the 
following: 

Bip = (Ns + N3)di3c(7i-73) (46) 

= N2N3[di3c(7i-73)]^ (4?) 

Since the maximum Bp for a 
given Bip occurs when = N 3 ^ 
equations (46) and (47) can be 
combined to give 



(48) 


The maximum Bp is plotted as the solid line on the parameter plane in 
figure 7- Since the minimum Bp is positive^ the system is, as previously 
indicated, stable. 

The stability of the system with errors in mechanization is determined 
next. Since a closed-form solution for the maximization of Bp does not 
exist, the equations were programmed on a computer. The maximum Bp was cal- 
culated for a fixed Bip over the range of parameters considered; that is, 

- 60 ^ ^ 7i^7s ^ +6o^, 1.11x10"^ ^ ^ 1 (equivalent to initial attitude 

errors as large as 5^), and mechanization errors of 0 ^ IA 7 i|, 1 A 73 | ^ 5*^* 

For gimbal angle errors of |A7 i| = |A7s| = 5^, figure 7 indicates an 
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insignificant difference in maximuin Bj) due to imperfections. The minimum 
Bj) is made positive by selecting the magnitude of the restricted region 
Since following inequality is derived from (43) 

Ibl < icot(7i-73)| . (49) 

where the minimum value of the cotangent is determined by the magnitude of the 
restricted region (A 7 j^). If |A7 i1 = -|A 73 |^ (49) can be written 

If [A/j^I = 10°, the right side of (50) is evaluated at I 71 - 73 I = 90° ±10°. 

The maximum gimbal angle error is therefore 

< 5 - 05 ° 

Therefore, B-n is always positive and the system is stable for any 
mechanization error less than 3 * 

Performance . - The performance of the system could be severely affected if 
large errors in mechanization existed and were not considered. For example^ 
assume that the system with perfect mechanization is designed so that the 
steady-state requirements are just met at the minimum loop gain^ which occurs 
on the boundary of the restricted region. If errors in mechanization actually 
do exists the system could be operating in the restricted zone and would not 
exhibit the required steady-state performance. Therefore^ the region in which 
gimbal angles cannot be commanded must be expanded by an amount equivalent to 
the maximum mechanization errors; thus^ the region of restricted operation is 

90° - (|A7r| + |8|^^) s |7i-7s| * 90° + (|A7 k1 + (51) 


Constant Processor With Tracker Pairs Mounted on 
Perpendicular Faces 

Constraints on the constants . - The stability of the system deriving its 
error signals from the constant processor for tracker pair 1,3 is determined 
in this section. The characteristic equation (CE) is 

CE = DET(i + GM^N) = 0 (52) 


where 
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rll 

1 

ri3 


0 

^S3 

^31 

0 

I’33 


or^ equivalently^ 

CE = (1 + sNi)(L + C^S + Cj^g^) = 0 (53) 

where 

= N2(r2iS7i - V 23 CJ 3 ) + WsCrsicyi + rsss/s) (54) 

Cj) = W2N3(r2ir33 - r3ir2s) 0 ( 71 - 73 ) (55) 

and the r. . are constants to be determined. 

-4 J 

Except for the definitions of the constants, equations (53) and (40) are 
the same. Therefore, the following two conclusions can be extracted from the 
analysis for the partial processor (imperfect mechanization). First, the ele- 
ments rii and ri3 can be equated to zero because they influence only the 
coupling between the control channels and not the stability of the system. 
Second, the stable ranges of variation of the parameters Cip and are the 
same as those of and Bq because the systems have the same parameter plane 

diagram. The diagram for the constant processor given in figure 8 shows not 
only the stability boundary illustrated in figure 7 sketch (e), but also 

the curve for ^ = 0.2. 

The parameters C,p and must, at least, be positive for the system to 
be stable. This requirement can be met only if the signs of the constants 
are controlled as a function of the outer gimbal angles 7^ and 73. Since 
N2 and N3 are positive, the following inequalities, which are the 22 and 33 
elements of the product matrix (see eq. ( 8 )), 


(rais7i - r23C7s) >. 0 

( 56 ) 

(r3iC7i + 1 - 33373 ) > 0 

(57) 


are sufficient to insure that C,ji remain positive. Since mechanization sim- 
plicity is desired, arbitrarily let r^i = 0. Inequality (56) is satisfied 
3^23 is negative since the gimbal range is -60^ 7i^73 ^ +6o^. It fol- 
lows from equation (55) that the parameter will always be positive if 

the sign of the constant rsx is controlled as follows. 

rsi ^ 0 if I71-73I $ 90° (58) 
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Figure 8. - Parameter plane - constant processor 
trackers 1^3* 




Because can be negative^ 

the following two conditions must 
be imposed on tqq to satisfy 
the inequality (57)* Firsts the 
sign of the constant r^s must 
be controlled as follows: 

^33 ^0 if 73 ^ 0 (59) 

or rs3S73 > 0 . Second^ the 
magnitude of Vss must be chosen 
so that the inequality of (60) is 
satisfied. 


kssl > 


I —I 

I ^^sjmax 

\ 3^31 



rsi 


(60) 


The region in which (6o) must be 
evaluated is shaded in 
sketch (f); that is^ the region 
in which r3i < 0. Also indi- 
cated on the sketch is the sign 
of the constants as a function 
of the gimbal angles (71^ 7s) 
and the region of restricted 
operation. Evaluating the 
inequality of (60) for the worse 
condition (71 = + 40 ° and 
73 = -60^) results in 

IrssI > 0.884lrsil (6l) 

Therefore^ the parameters Cj^ 
and Grp are positive if = 0^ 

rs3 < 0 ^ and inequalities ( 58 )^ 
(59)^ and (61) are satisfied. 

The parameters Cp and Cip 
are positive but the system is 
not necessarily stable. The con- 
stants (r^s^ ^31? ^33) 3 ms t be 
chosen so that the range of Cj^ 
and Cp is contained within the 
stable operating region indi- 
cated on the parameter plane 
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(fig. 8). Searching for a set of constants over an unrestricted range that 
will provide a stable system is impractical. The range of constants over 
which the search must be made is reduced considerably by investigating only 
those sets that provide the desired steady-state performance. 

Steady- state performance . - The equations describing the steady-state 
pointing error and the steady -state drift rate error for the system using the 
constant processor are derived in appendix C and are given by 


v^ss ^ss 



ttss^so) + (-^32^20 + ^22^30) ~ (62) 


2 / 2 

^ + h22 ^ DR (63) 


where 

H = h^j = for tracker pair 1;,3 

Hp = (r^irss - r23^3i) c(7i“73) 
tj_ = external torque about the ith axis 


^10 


n j 

^0 “ ^lO 



• J 

^20 


^0 " ^ ^W20 

kso 


GO 

^0 W30 


hpQ is proportional to the initial momentum about the ith axis. If 
Hj) and the elements of the H matrix are substituted for trackers 1^3 
(eq. (8)) and = Oj equations (62) and (63) can be written as quadratic 

equations in rps as follows: 


{h|o - [(PE)Krsic(7i-73) ]^}r|3 - {2[r33 - r3is(7i-73) ]h2oh3o}3T23 


+ [^ii + ^33 " 2r3ir33s(7x“73) ]hfo ^ 0 (64) 


{^ 3 '^ - [ (DR)IiKr3ic(7i-73) ]^3r|s - {2[rs3 - r3is(7i-73) ]t2ts]r23 


+ [r|i + r %3 - 2r3ir33s(7i-73 ) ]t2" 


(65) 
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Sets of constants ^33) 

can then be calculated so that 
the system will satisfy each of 
the individual steady-state 
req.uireinents for the most severe 
conditions. These conditions are 
listed in table III. The accept- 
able range of rss for a given 
r3i and rss is plotted in fig- 
ure 9- The solid lines of con- 
stant |rsi| indicate the minimum 
1^23 1 Tor a given Irssj to sat- 
isfy the drift rate specif ication. 
The dashed lines provide the same 
information for the pointing accu- 
racy specifications. For example^ 
if 1^31 1 = 2.0 and Irsaj = 3.O, 
the drift rate accuracy requires 
that jr^sl ^3.65 and the point- 
ing accuracy requires that 
1^2sl ^ 3»32. In general, the 
jr^sj that satisfies the drift 
rate specification also satisfies 
the pointing accuracy 
specification. 

Stability . - The stability of 
the system is now investigated by 
calculating the range of the vari- 
able parameters (C^, Cp) for a 
specific set of constants chosen 
from the acceptable region. The 
acceptable region is defined as 
the region that satisfies not 

only the steady-state conditions but also the inequality of (6l). The hatched 
boundary in figure 9 is plotted for the condition that 1^33 j = 0.884 |r3ij. 
Therefore, for a given rss, acceptable sets of constants are chosen from the 
region to the right of the hatched boundary and above the constant bound- 

ary. The sets of constants investigated are listed in table IV(a). Sets 1, 

2, and 3 provide a stable system whereas sets 4 and 5 result in an unstable 
system. The "worst case" results of the stability calculations for sets 1 and 
2 are plotted on the parameter plane in figure 8. The sets are also shown on 
figure 9 with an indication of system stability (the A indicates a stable 
set and the "o" indicates an unstable set) . The figure shows that sets of 
constants that assure stable operation begin to occur when the left side of 
(57) becomes significantly positive. 



Figure 9. - Steady-state gain requirements - 
constant processor - trackers 1,3- 
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Constant Processor With Tracker Pair Mounted 
on Parallel Faces 


Constraints on the constant s, - The stability and performance of a system 
using trackers 1 and 2 is investigated next. The characteristic eq.uation is 
identical to equation (52) except that and W are defined for trackers 1^2 
(table II). Therefore^ the form of equation ( 53 ) trackers 1^2 is 


CE = (1 + gNi)(l + E^g + Ej^g^) = 0 (66) 

where 

E,j, = N2(q2iS7i - 123^72) + %(lsiC7i - Q.33<^72) (^7) 

E-q = W2Ns ( 923131 ~ I21I33) 2(71-72) ( 68 ) 



Figure 10.- Parameter plane - constant processor - 
trackers 1,2. 


As in the previous case^ the con- 
stants qii and qxs are equated 
to zero. Also^ the parameter 
plane diagram for the system^ 
using the constant processor for 
trackers 1,2 (E,p, ^ identi- 

cal to the diagram for the imper- 
fect mechanization of the partial 
processor (sketch (e) or fig. 7) 
and is given by figure 10. 

A necessary condition for 
stability is that E,j,, > 0. 

The parameter E,j, is positive if 
the following inequalities are 
satisfied. 

q^is/i - > 0 (69) 

<l3iC7i - «l 33 C 72 > 0 (70) 

To simplify the mechanization, let 

l2i = l23- inequality of (69) 

is satisfied if the constants q^i 

and q23 are controlled as 
follows. 

123 = I21 if (71-72) ^ 8 

(71) 
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The parameter Ep is then positive if 


(Q31-Q33) > 0 (72) 

It remains to insure that the inequality of (70) is satisfied. It follows 
that (70) and (72) are satisfied if qsi and q^s are chosen as indicated 
below. 


<331 > <l 33 > 0 

if 

l 7 i 

1 < l 7 a| 

(73) 

<lS 3 < <331 < 0 

if 

l 7 i| 

1 > I72I 

( 7 ^) 

are then positive. 

The 

sign of the oonstants 

as well 


raraine uera tuiu. jliji 

as the region of restricted operation is plotted as a function of the outer 

gimbal angles in sketch (g)j also 
031^^3 <0 is shaded area and q^x >0 is 
cross hatched area. 


Steady-state per form ance . - The equations 
for steady-state pointing accuracy and steady- 
state drift rate accuracy are given, respec- 
tively, by equations (62) and (63) where H 
^2 and Hj) are now defined for trackers 1 and 2; 
that is 


H = h^j = for tracker pair 1,2 

= q2i(q3i~<l33)s(7i-7a) 

Substituting the elements of the H matrix 
into equations (62) and (63) yields 



Sketch (g) 


{ [ PEK( qsi -qs 3 ) s (/i -72) ] - 2[1 - 0(71-72) ]h 3 o}q 2 i 

+ ( {2(qsi+q33) [1 - ^( 71 - 72 ) ] ]h2ol^3o)^i ■ t ^ii + “ 2q3iq33c(7i-72) ]h|o ^0 


( 75 ) 


[DRK:(q 3 i-q. 33 ) 3(71-72)]^ - 2[i - 0(71-72)] (^) [ q|i 


+ ( {2(q3i+q33) [1 “ 0(71-72) ]} t” ^ j 




t 2 t 3 

r 

2 


- [^31 + 0.33 - 2q3iq33c(7i-72) ] (^ 7 ^ ^ ^ ( 76 ) 
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(These equations are also given by equations (C23) and (C24) in appendix C.) 
The acceptable region of q^i for a given qsi and qss can then be deter- 
mined for the most severe condition from equations (75) and (T^)- If the 
equations are evaluated at the conditions listed in table the most severe 

condition is determined by scanning the range of gimbal angles and by control- 
ling the sign of the input parameters (tj_ and hj^o)* constant (q^i) that 

satisfies the most severe condition is plotted in figure 11^ where the solid 




^31 

Figure Li.- Steady- state gain requirements - constant processor - trackers 1,2. 

lines of constant qss indicate the minimum q^i necessary to satisfy the 
drift rate specification. The dotted lines of constant q_ss indicate the 
minimum required to satisfy the pointing specification. The acceptable 

q^i for a given qsi lies above the constant qss curve- The diagram shows 
that the sets of constants which satisfy the drift rate specification also 
satisfy the pointing accuracy specification- The symbols A and o in fig- 
ure 11 represent sets of constants that satisfy the steady-state requirements 
and will be discussed in more detail in the following section. 

Stability . - Stability of the system is determined for sets of constants 
that satisfy the steady-state requirements. System stability was investigated 
for the sets of constants given in table IV(b). As indicated in the table^ 
sets 1, 2, and 3 provide stable operation of the system whereas sets 4 and 5 
provide an unstable system over the range of parameters investigated 
( -60° ^ 71^72 ^ +60^ and 1.11x10'“^ ^ ^ 1) • The degree of stability for 
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sets 1 and 2 is indicated by the parameter plane diagram of figure 10 where 
the plotted points represent the solution for the worst case. The sets of 
constants are also shown on figure 11 with an indication of the system stabil- 
ity. The A and o represent^ respectively^ a stable set and an unstable set. 
The superscript indicates the particular set of constants. 


SIMULATION 


The attitude control system was simulated on the analog computer to eval- 
uate the performance of the system using each of the processors. The charac- 
teristics of the satellite and the control system are given in appendix D. 
Also, the star tracker gimbal angles and the characteristics for each 
processor are included in the appendix. 

The purpose of the simulation was twofold. First, it was desired to val- 
idate the results of the stability analysis obtained from the digital com- 
puter. Second, it was desired to compare the performance of the system with 
the different processors. The simulation included the saturation type non- 
linearity introduced by the motor but excluded the gyroscopic coupling torques 
due to motor rotation. The geometry transformation was described by the 
linearized equations (N matrix) . 


Ideal Processor 

The system using the ideal processor to derive the attitude control sig- 
nals was simulated on the analog computer. The system was designed so that 
the linear damping and natural frequency are, respectively, ^ = 0.7 and 
copj ^ 1.2. From the parameter plane (fig- 5)^ the variable parameters are as 
follows: 


B1 = — = 2.0 
^2 

B2 = - 14.32 where (N^ = l) 

As is to be expected from the relation in figure 6, the system then provides 
the desired steady-state performance for external torques less than 
4750 dyne -cm. 

The transient response of the system with the ideal processor is given in 
figure 12. The attitude error about each control channel is plotted twice so 
that the system response can be observed not only for large deviations but 
also for small deviations in the region near equilibrium. The figure shows 
that the control channels give identical responses to initial errors in 
position and are decoupled. 
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Figure 12.- Transient response of system using the ideal processor 


Partial Processor 


The system using the partial processor vas designed so that it exhibited 
the desired performance over the ooraplete range of operating conditions. The 
minimum linear loop gain is chosen so that the desired steady-state 

performance is achieved for external torques less than or equal to 

2000 dyne-cm (fig. 6 ); thus^ = 5* 96- From equation (37)^ B2max = 3^-^ 

for l 7 i- 7 sl = 80*^ which corresponds to a restricted region (A 7 j^) of 10^. The 
compensator time constant (b 1 = is chosen so that satisfactory perfor- 
mance is achieved over the complete range of B2. If B1 = 1.25, the damping 

ratio of the system varies from ^ = 0.25 to ^ ^ 0-70 as observed from the 

parameter plane (fig. 5 )* 

The transient response of the system is indicated in appendix D (fig- 
ures 13 , l4, and 15)* The response of the system to initial attitude errors 
of 5^ about all control axes is shown in figure I 3 when the loop gain is mini- 
mum (B2 = 5 . 96 ) and in figure 1^ when the loop gain is maximum (B2 = 34.4). 

The difference between the two response curves is negligible because the lin- 
ear range of the saturation type nonlinearity is so small that it functions 
like a relay; thus, the system functions in the nonlinear region most of the 
time. 


The coupling between control channels is observed by comparing the three 
transient response runs displayed in figure 15* Each run shows the response 
of the system to an initial attitude error of 5 ^ about a single control axis. 
Comparison of the runs indicates that pitch and yaw motion is coupled into 
the roll channel. Also, the pitch and yaw channels are observed to be 
independent. 


Constant Processor 

The system using the constant processor for tracker pair 1,3 and tracker 
pair 1,2 was also simulated. For a given set of constants^ the gimbal angles 
were chosen in each case so that the system was operating under the most 
severe conditions . 

The transient response of the system using the constant processor for 
tracker pair 1,3 is given in appendix D (figs. 16 and 17)* The constants for 
the processor are r^s = -4.25, l^si| = 2.0, and \tss\ = 3* 5* Figure 16 shows 
runs in which an initial attitude error is imposed about a single control 
axis. The run on the left of the figure shows that roll motion is not coupled 
into pitch and yaw whereas the remaining two runs show that pitch and yaw 
motions are coupled into the other two channels. Figure I 7 shows the response 
of the system when an initial attitude error is imposed simultaneously on all 
three control channels. 

The response of the system using the constant processor for tracker pair 
1^2 is given in appendix B (figs. I 8 and 19)* The response in both figures 
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is for the same set of constants (jq^il = Q31 = 


+ 7-0 

- 6 . 0 ’ 


and qs3 


+6.0v 

- 7 - 0 ^ 


but for different gimbal angle values. Figure 18 shows the behavior of the 
system for gimbal angle values that provide the most unsatisfactory response. 
The yaw channel response clearly indicates that the coupling term dominates 
until the attitude error from the pitch channel sufficiently decreases. Fig- 
ure 19 shows the behavior of the system for gimbal angle values that provide 
significantly improved response. 


CONCLUSIONS 


Three methods of processing star tracker information to derive attitude 
control signals have been developed. They are: (l) ideal processor, ( 2 ) par- 

tial processor, and (3) constant processor. The system using each of the pro- 
cessors exhibits acceptable transient and steady-state performance. Also, 
each of the processors is simple enough to be mechanized on board current 
spacecraft. 

The star trackers are used in pairs to estimate the attitude errors of 
the satellite via the processors. The simplest form of the processors is 
obtained when the paired trackers are mounted to the vehicle so that their 
outer gimbal ax:es are parallel and when two inner and one outer gimbal error 
signals are used to derive the attitude control signals. The stability of 
the system is then independent of the inner gimbal angles of the star trackers. 

An easily applied criterion for selecting stars (by ground operation) to 
avoid the indeterminant condition exists if the star trackers are used as 
suggested. The indeterminant condition implies that the attitude of the vehi- 
cle is unobservable; that is, the gimbal error signals being used to estimate 
the attitude errors are dependent. The condition occurs when the plane formed 
by the optical axes of the paired trackers contains the outer gimbal axes of 
the trackers. A region of restricted operation is established in the neigh- 
borhood of the indeterminant condition to insure that desired performance is 
maintained. 

The ideal processor is the most complex of the three processors from the 
point of view of implementation because such trigonometric functions as the 
sine and cosine of the difference of gimbal angles and the tangent of gimbal 
angles must be mechanized. However, the control signals derived from this 
processor are independent and therefore provide a well-behaved system. The 
desired system performance is achieved over a commanded gimbal range that is 
limited only by the physical nature of the star tracker. 

The partial processor is much simpler to mechanize than the ideal pro- 
cessor. Only the simple trigonometric functions such as the sine and cosine 
are mechanized with resolvers. The processor provides two control signals 
that are independent. Also, the required performance of the system is 
achieved over a commanded gimbal range that is again limited only by the 
physical nature of the star tracker. 
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Errors in mechanization of the partial processor have a negligible effect 
on the stability and performance of the system. To avoid problems due to 
errors in mechanization, it is only necessary to insure that the errors do not 
allow the system to operate in the region of restricted operation. Therefore, 
the restricted region must be increased by an amount equivalent to the error 
in mechanization. 

The constant processor is the simplest of the three processors. Only 
amplifiers, potentiometers, and relays are required to mechanize terms that 
are constant in magnitude but vary in sign as a function of the outer gimbal 
angles. Performance of the system using the constant processor is acceptable 
over a gimbal angle range of ± 60 ^. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif. , 9^035, Nov. 29, 196? 

125-19-03-06-00-21 
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APPEKDIX A 


DERIVATION OF EQUATIONS RELATING THE MOTION OF 
THE TRACKERS AND VEHICLE 


The relationship between the motion of the trackers and the motion of the 
vehicle is derived in this appendix- The relationship is obtained by express- 
ing the motions as angular velocities and exploiting the simple transformation 
properties of angular velocity vectors. 

Assuming that the tracker tracks its star without error^ then the angular 
velocity of the vehicle in inertial space is the same as its velocity relative 
to the trackers. The inertial velocity can be expressed in terms of compo- 
nents along orthonormal axes fixed in the vehicle as coiyl^ + 

The parameters of the tracker are needed to express the velocity of the 
vehicle relative to the tracker. 

The sketch on the right of figure 2 shows the angles associated with 
tracker 1. These angles^ with appropriate subscripts^ are also associated 
with the other trackers. Each angle a relates a "line of sight" coordinate 
system {l) = (l^^ 2^^ 3^) to a coordinate system (t) = (l^> 2^^ 3^) fixed to 
the tracker in such a way that 1^ == ft. The (t) coordinates are related to 
a (y) = (l^^ 2^ y 3^) coordinate system by p so that 2^ = 2^. Similarly^ (y) 
is related to an (x) = (l^> 2^^ 3^) coordinate system by 7 so that 3^ = 3^* 
The angles are positive when generated by right-hand rotations about their 
generating axes. Whereas the (z), (t)_, and (y) coordinates and their angles 
P^ and 7 move as the vehicle moves, the (x) coordinate is fixed to the 
vehicle in a way that varies from tracker to tracker. The (x) coordinate of 
each of the four trackers amounts to a permutation of the vehicle’s coordinate 
system (v): thus (x) = P(v), where 



The inertial angular velocity of the vehicle in (x) coordinates is given 
as + ^3x3^* The fundamental relationship which states that 

this velocity is equal to the velocity relative to the tracker can be 
expressed in the form 


“ixl^ + “2x2^ + Wsx3^ = + P2^ + 73^ 


(A2) 


The following equations are obtained from the definitions of the coordinate 
systems: 



(A4) 


y t 

The 2 equation and the 1 equation are obtained by taking the inverse of 
equations (A3) and (A4) and are substituted into equation (A2) to give 



(A5) 


Inverting this 


transformation 


to solve for 


the tracker parameters gives 



(a6) 


Equation (a 6) holds for all the trackers^ which differ from each other 
by the way they are mounted on the vehicle. The manner of mounting is 
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expressed by the permutation matrices in equation (Ai). Postmultiplying the 
transformation matrix in equation (a 6) by for example^ relates the num- 

ber 1 tracker parameters to the components of inertial angular velocity in 
vehicle coordinates : 



C7i 

S7icPi 

-C7ispi 



(AT) 


The rate transformation (R) for each of the trackers in figure 2 is obtained 
in a similar manner and is given in table I. 


4l 



APPEKDIX B 


DERIVATION OF THE CHARACTERISTIC EQUATION FOR THE SYSTEM USING 


THE PARTIAL PROCESSOR WITH IMPERFECT MECHANIZATION 


The equations describing the effect of errors in mechanization of the 
partial processor on the stability of the system are derived in this appendix. 
The errors are assumed small enough so that only the first-order terms are 
significant. 

The partial processor for tracker pair 1,3 "when perfectly mechanized is 


M-c 


0 10 
^13^73 0 “di3C7i 

^^13^73 0 di3S7i 


(Bl) 


It is assumed that the constant terms are perfectly mechanized and that errors 
exist only in the gimbal angles 7 i and 73 . The errors result from either a 
misalinement of the resolvers used to mechanize the trigonometric functions or 
an inaccurate measurement of the gimbal angle. Gimbal angles with errors are 
represented by + ^ 7 i* partial processor with errors in mechanization 

is 


M- 


PE 


di3s(73+A7s) 

^■ 132 ( 73 +^ 73 ) 


1 0 
0 -di3c(7i+A7i) 

0 d-i3s(7i+A7i) 


The characteristic equation is 


CE = DET(l+GMpgN) = 0 


where 


G = 


Si 0 0 

0 g2 0 
0 0 g 3 


Si = 


WiK(TiS + 1) 
s(t2S + 1)(Tj„S + 1) 


(B2) 


(B3) 


(b4) 
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I 


Define the following 


so that 


g 0 0 


G = 0 g 0 , g 

0 0 g 


K(tis + 1) 

s(t2S + 1)(t„s + 1) 


B - DpMpgN - 


Ni 0 0 


0 No 0 


0 0 Ns 


hii hia hi3 

0 bga t>23 

0 t>32 t)33 


G = GD. 


The characteristic equation Is 


CE = DET(I + GB) = 0 


CE = (1 + gNi)(l + B,jg + Bpg^) = 0 


where 


Bip = ^22 B33 


(BIO) 


Bj-) = 1*22^33 " 1*23^32 


(Bll) 


The B matrix for the partial processor for tracker pair 1^3 Is 


-NiC/itPi 


NiS7it3i 


B = 0 {N2<ii3[s(73+A73)s7i + c(7i+A7i)c73]} {N2<ii3[s(73+A73)c7i - c(7i+A7i)s73]} 

0 {N3<ii3[c(73+A73)s7i - s(7i+A7i)c73]} {N3<1 i3[c(73+A73)c7i + s(7i+A7i)s73]} 

(B12) 
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If A 7 j is small enough so that the approximations cA7p = 1 and sA7^ = AXj^ 
are valid, equation (B12) can he written as 


Ni 


-NiC/itPi 


NiS/itPi 


B = 


0 {W2di3[c(7i-73) + s7iC73(A7s-A7i)]} {NgdisE ( 073071 )A7s + (s7iS7s)A7x]} 


0 


(-N3di3[(s73S7i)A73 + (c7iC73)A7i]} 


{Nad isC 0 ( 71 - 73 ) 


- C7iS73(A73-A7i)]} 


(BI 3 ) 


By using equation (B13)^ the parameter Bij (BIO) can now he written as 


B^ = (W 2 + Ns)di3c(7i-73) + di3(N2S7iC73 - W 3 C 71 S 73 ) (A 73 -A 71 ) (B14) 

Likewise, the parameter B^ (B11) can he written as follows. 

Bj) = W 2 N 3 df 3 {[ 0(71-73)]^ + 6 [c( 7 i- 73 )s( 7 i- 73 )] - 5 ^[ S73C71C73S71] 

+ [(c73C7i)A73 + (s7iS73)A7i][(s73S7i)A73 + ( C73C7i)A7x] } (B15) 

where 


6 = (A73-A71) 

Equation (B15) reduces to the following when only the first-order error terms 
are considered significant. 

Bp = N2Wsdf3c(7i-73 ){c( 7 i- 73) + S[ 3(71-73 )]} (B16) 

Equations (B14) and (B16) define the parameters that are used in the charac- 
teristic equation (B9) to determine the effect of errors in mechanization on 
the stability of the system. 
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Sketch (h) can be reduced to sketch (i) 



where 



(C3) 


Tm J 




(C4) 

(C5) 


The steady-state performance is determined by writing the output parameters 
as a function of the input parameters and using the final value theorem; that 
is^ 



The steady-state errors are determined as follows. The steady-state 
position error caused by an initial angular velocity of the vehicle and an 
Initial angular velocity of the motor is 
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Ej^ = determinant of H 

The velocity error for a constant external torque about all three axes 
[Ti(s) = ti/s] Is 



The steady-state pointing error and drift rate error are defined^ 
respectively^ by equations (C9) and (CIO). 

(Cll) 

(C12) 


Pointing error = (0^^ + ^ PE 

Drift rate error = (^1 + ^ DR 
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where PE and DR are, respectively, the maximum allowable pointing error and 
the maximum allowable drift rate error. Use of equations (C7) and (CIO) 
allows equations (Cll) and (C12) to be written as follows; 



[ (hsshso 


has^So)^ + (-^32^20 


^ PE 


(C13) 







h22 ^ 


2 Hi/ a 


h DR 


(Cl4) 


The system will perform as required in steady state if the inequalities of 
(CI 3 ) and (Cl4) are satisfied. The equations are general since the constant 
processor for a particular tracker pair defines the H matrix. The steady- 
state performance of the system using the constant processor is now determined 
for the following two specific cases. They are: (l) tracker pair 1^3 

( 2 ) tracker pair 1^2. 


Tracker Pair 1^3 

The H matrix for tracker pair 1^3 is 


H = M^N 


1 (niiS7i-c7itPi-risC73) (rnC7i+S7itPi+ri3S73) 

0 (r2iS7i-r23C7s) (i’2iC7i+n2sS73) 


0 (r3iS7i-r33C7s) 


(n3iC7i+r33S73) 


(C15) 


and 


= (r2ir33-r23r3i)c(7i-73) 


(CI 6 ) 


With T 2 i = 0 and equations (C15) and (Cl6), the pointing error equation (C13) 
becomes 


J [^31+^33~'^^3X^33^{7i.-73) ]hfo “ { 2^23 [ ^33 -^SlS ( 7l "73 ) 1 l^^aO^SO 


^l3hio 


K[ -r2sn3ic(7i-73)] 


^ PE 
(CIT) 


and must be evaluated at the most severe conditions. If it is assumed that 
= 10*^ and \h^o\ = j^soj^ (d?) is evaluated at the following conditions. 
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^21 = 0 


7 i -/3 = + 80 ° 


^20 ~ “^30 ^31 > 0 (Cl8) 

r23 < 0 rs3 < 0 


In order to determine acceptable sets of constants, (CI 7 ) is written as a 
quadratic in r^s as indicated by 

{bfo - [ (PE)Kr3ic(7i-73) - {2[r33-r3is(7i-73) ]h2ob.3o}^23 


+ [rfi + r |3 - 2rsir33s(7i-73) ]h|o ^ 0 (C19) 

and is evaluated subject to the conditions listed in (Cl 8 ). The allowable 
range of r 23 for a given rsi and Tss that will satisfy the pointing accu- 
racy specifications over the complete range of operating conditions can be 
determined. 

The constants required to satisfy the drift rate specifications are 
determined in a like manner. In fact^ from (C13) and (Cl4)^ the corresponding 
drift rate inequality can be written from (C19) and is 

{t3^ - [(DR)KIir3ic(7i-73)]^}r|3 - {2[r33-r3is(7i-73) ]t2t3]r23 

+ [rfi + r|3 - 2r3ir33s(7i-73) ]t 2 ^ ^ 0 (C20) 

Equation (C20) is evaluated subject to the conditions listed in (Cl 8 ) and 
ta = -^ 3 . The acceptable range of ras for a given and that will 

satisfy the steady-state drift rate error for all operating conditions can 
then be calculated. 


Tracker Pair 1^2 

The H matrix for tracker pair 1^2 is 


H = 


1 (tliiS7i-c7it|3i-qi3S72) 
0 021 ( 371 - 372 ) 

0 ( 031371 - 0 . 33372 ) 


(0iiC7i+S7it3i-qi3C72) 

021 ( 071 - 072 ) 

(O 31 C 71 -O 33 C 72 ) 


(C21) 


where 
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l2i = 'las 


He = qai(qsi-'ls3)s(7i-7a) 

The elements of the matrix in equation (C21) can be substituted into 
equation (C13) to provide the following equation 

([ill + lis " 2qsiqssc(7i-72) ]hfo - [SqaiCqsi+lss) [1 - c (7i -7a) 3 iHaoH-so 

+ {2q|i[l - c(7i-72)]}h|o) ' [Kq^iCqai-qaa) s ( 71 - 72 ) ] ^ ^ PE (C22) 

The inequality of (C22) is written as a quadratic in (C23) so that sets 

of constants can he determined that will satisfy the steady-state pointing 
accuracy. 

[[ (pE)K(q 3 i-q 3 s) 3(71-72) - 2[1 - c( 7 i -72 ) Jhloloji 

+( (2(131+133 ) [1 " ^(71-72) ] }h 2 o^ 3 o) I21 

- [ill + ^3 ” 2q3iq33c(7i-72) ]hfo ^ 0 (C23) 


The most restricted range of 121 for a given 131 and 133 is determined if 
(C23) is evaluated at the most severe conditions. 

As before^ the inequality for drift rate error can be written directly 
from (C23) and is 


[ (dR)k( 131-133)3(71-72) - 2[1 - 0(71-72)] (^) 


+ ^{2(q3i+q3s) [1 - c(7i-7a)]} ^ lai 


- [ill + I33 - 2qsiq33c(7i-72) ] ^ ® (C2l) 


Again the parameters are chosen so that the worse case is obtained; that is, 
it is necessary to choose the sign of t2 and t3 and the magnitude of (71-72) 
that results in the most severe condition. The results of the calculations 
are plotted in figure 11. 
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APPEMDIX D 


CHAEACTERISTICS OF IHE SATELLITE 


Vehicle characteristics 

Ii = = 1952 kg-m® (l 440 slug-ft^) i 

"ij = ° 

Motor characteristics 

Stall torque = Tg = 0.0353 N-m (5 in. -oz) 
Time constant = = 76.8 sec 

= 0.1041 N-m-sec (O.O768 ft-lh-sec) 
Compensator characteristics 
Ideal processor 

= 268^000 


= 1. 2, 3 


c radians 

Ti = 5-27 sec 

T2 = 0.527 sec 

Partial processor 

K = 268,000 — ^ 

c ^ radians 

Ti - 8.0 sec 

t 2 = 0.8 sec 

Constant processor 

Kc - 268,000 

Tj_ = 5*0 sec 
T2 = 0.5 sec 
Processor and geometry 
Ideal processor 

Unity feedback 
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III I III I Mill I II I lllMaiB* 


Partial processor 

jdxsi = 2.4 
Pi = 30° 

7i = 60° 

Figure 13 

71-73 = 80° "where 73 = -20° 
Figures l4 and 15 

7i-7s = 0° where 73 = 6o° 
Constant processor 

Tracker pair 1,3 - figures 16 and I7 



-4.25 

^31 

= 2.0 

I’33 

= -3-5 


7i = 

-35° 

Pi 

= 30° 

73 

O 

LT\ 

-it 
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Tracker pair 1^2 
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Figure l8 
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Figure 15»- Transient response of system using the partial processor - initial position error 

about single axis to show coupling. 
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Figure 19* - Transient response of system using 
the constant processor for trackers 1,2 - 
arbitrary gimbal angles. 
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TABLE I. - TRANSFORMATION R RELATING THE STAR TRACKER 


GIMBAL RATES AND THE VEHICLE INERTIAL RATES 
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TABLE II.- TRAWSFOEMATIONS DESCRIBING THE GEOMETRY (N) AND THE PROCESSORS 

(M, Mp^ M^) FOR PAIRS OF TRACKERS 
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TABLE III. - DATA REQUIRBID FOR EVALUATING 
STEADY -STATE PERFORMANCE OF SYSTEM 
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TABLE IV. - SETS OF CONSTANTS FOR THE CONSTANT PROCESSOR 
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